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Abstract In this study we use a numerical simulation of an artificial coronal 
mass ejection (CME) to validate a method for calculating propagation direc- 
tions and kinematical profiles of interplanetary CMEs (ICMEs). In this method 
observations from heliosphcric images are constrained with in-situ plasma and 
field data at I AU. These data are used to convert measured ICME elongations 
into distance by applying the Harmonic Mean approach that assumes a spherical 
shape of the ICME front. We use synthetic white- light images, similar as observed 
by STEREO-A/HI, for three different separation angles between remote and in- 
situ spacecraft, of 30°, 60°, and 90°. To validate the results of the method they 
are compared to the apex speed profile of the modeled ICME, as obtained from a 
top view. This profile refiects the "true" apex kinematics since it is not affected 
by scattering or projection effects. In this way it is possible to determine the 
accuracy of the method for revealing ICME propagation directions and kinemat- 
ics. We found that the direction obtained by the constrained Harmonic Mean 
method is not very sensitive to the separation angle (30° sep: (j) = W7; 60° sep: 
= W12; 90° sep: cj) = W15; true dir.: EO/WO). For all three cases the derived 
kinematics are in a relatively good agreement with the real kinematics. The best 
consistency is obtained for the 30° case, while with growing separation angle the 
ICME speed at 1 AU is increasingly overestimated (30° sep: AVarr ~ —50 km 
s-\ 60° sep: AKn- ~ +75 km s-\ 90° sep: AKrr ~ +125 km s-^). Especially 
for future L4/L5 missions the 60° separation case is highly interesting in order 
to improve space weather forecasts. 
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1. Introduction 



Coronal Mass Ejections (CMEs) are the most powerful eruptions on the Sun 
and are currently observed by a relatively large number of instruments. Space- 
based coronagraphs are located at three different positions to monitor the solar 
corona from aU sides. SOlar and Heliospheric Observatory (SOHO) LASCO 
coronagraphs (C2 and C3: [Brueckner et al. 1995) at Li provide white-light 
images up to a distance of 32 Rq. The two spacecraft of the Solar TErrestrial 
RElations Observatory (STEREO: Kaiser et al. 2008) move apart by 45° per 
year and are currently at a position of more than 120° eastern and western 
from Earth, respectively. As part of the Sun-Earth Connection Coronal and 



Heliospheric Investigation suite (SECCHI: Howard et al. 2008 Harrison et al. 
[20C)9| ) the CORl and C0R2 coronagraphs aboard the two spacecraft STEREO- 
A and STEREO-B provide continuous white- light observations up to 15 R©. 



Combining coronagraphs with the Heliospheric Imagers HIl and HI2 (Eyles 
et al.\ ,2009), which are also part of the SECCHI suite aboard STEREO, the 



entire propagation of an ICME from Sun to 1 AU is covered. Both His are off- 
pointed from the center of the Sun and cover a field of view (FoV) of 4° to 24^ 
and 18.7° to 88.7° elongation, respectively. 

Such wide FoVs make the interpretation of extended structures in these im- 
ages rather difficult and Thomson scattering effects fBillings 1966 Hundhausen 



1993) as well as effects due to line of sight (LoS) integration have to be taken 



into account. Methods were developed in order to investigate the propagation 
of CMEs beyond the typical FoV of coronagraph images, which are up to about 
4° (C0R2) or 9° (LASCO C3) elongation. These methods are based on specific 
geometrical assumptions about the shape of the ICME front. Some assume a 
point-like structure as the Fixed-0 method (F(/): Sheeley et al. 1999 Rouillard 
et al. |2008) or the triangulation method developed by Liu et al. (20101 that 



uses stereoscopic observations. The Harmonic Mean method (Lugazl 2010) and 



the Self Similar Expansion method (Davies et al. 2012) assume a circular front. 
Fixed-(/), Harmonic Mean (Lugaz 2010 Mostl et al. 20111 and the Self Similar 



Expansion method (Mostl and Davies 2012 1 are also used as fitting methods by 



assuming constant ICME speed and direction. In this way, these methods can 
be used to forecast the ICME arrival at 1 AU. An alternative forecasting tool 



is the drag based model, developed by Vrsnak et al. (2012). It is based on the 



idea that the aerodynamic drag in the ambient medium is responsible for the 
adjustment of the ICME speed to the solar wind speed. 

ICMEs are not only remotely observed but also measured in-situ, among 



others, from both STEREO spacecraft (IMPACT: Luhmann etal\ 2008), from 



Wind (Ogilvie et al. 1995) and the Advanced Composition Explorer (ACE: Stone 



et al. 



1998). ACE and Wind are both located at the Li point. In order to 



guard against misunderstandings, we use the term ICME as defined by |Rouillard 
(2011), i.e. an ICME includes the shock (if present), the dense sheath region and 
the magnetic flux rope. Furthermore we use this term for white- light observations 
beyond the coronagraphic FoV as well as for the in-situ structure. Using plasma 
and magnetic field data from in-situ observations, in addition to remote sensing 
of ICMEs, we are able to constrain the results derived from measurements of 
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white-light images ( |Mostl et aL]|2009a[[2010||Rollett eram2012[ ). In the following 
we refer to the method as "constrained" Harmonic Mean method (CHM). The 
kinematical profiles resulting from the CHM method may also be compared to 
the results of the drag based model and serve for studies of the interaction of 



ICMEs with solar wind strucures (Temmer et al. 2011[ 2012) 



MHD simulations of ICMEs evolving in interplanetary (IP) space have already 
been used to validate the YiKed-cj) fitting and Harmonic Mean fitting methods 
(Lugaz, Roussev, and Gombosi, 20111. By using synthetically derived ICME 
density images and time-elongation plots (Jmaps: Davies et al. 2009), measured 



from different viewing angles, we aim to validate the CHM method and present 
a detailed error analysis in order to improve the forecasting of Earth-directed 
ICMEs. We construct the Jmaps for different separation angles between the 
remote and the in-situ spacecraft observing the ICME, namely 30°, 60°, and 
90°. The 60° case is in particular interesting, as it simulates the observational 
capabilities for L4/L5 missions, which are currently under discussion to optimize 
space weather predictions ( jGopalswamy et al. 2011). 



2. MHD Simulations 



Simulations give us the opportunity to study the kinematical behavior of ICMEs 
from an ideal location as observer, namely top view along the Sun-Earth line, 
which enables us to directly measure the time-distance profile of the ICME 
evolution. This is in contrast to how we actually observe ICMEs, i.e. from 
a side view involving LoS and scattering effects, from which we measure an 
elongation angle that needs to be converted into radial distance by using different 
conversion methods. From the model we obtain both views, and use these two 
capabilities to test the efficiency of the CHM method to derive the kinematics 
of an ICME measured from synthetic Jmaps as observed by STEREO-A and 
to calculate ICME propagation directions. Synthetic Jmaps are produced from 
three different observing positions, i.e. 30°, 60°, and 90° between ICME apex 
and remote spacecraft in order to test the influence of the observing angle to the 
resulting directions and kinematics. 

The numerical modeling of the simulated ICME is done by using the Space 



Weather Modeling Framework (SWMF: |T6th et a^.[[2005| . Along the Sun-Earth 
line the grids of the model vary in spatial resolution from 0.5 R© close to the Sun 
to 0.25 Rq in the middle and 0.125 R© close to Earth. Synthetic in-situ data, i.e. 
magnetic field, plasma density and plasma bulk velocity, are produced for the 
location of Earth. This is done using the semiempirical MHD model developed 



by Cohen et al. (2007), which applies the empirical Wang-Sheeley-Arge model 



(WSA: [Wang and Sheeley[ |1990[ |Arge and Pizzo[|2000| and a flux ro pe CME 
initiated at the solar surface (similar setup as in Lugaz et al. 2009b). During 



the last 16 hours before the shock hits Earth the temporal resolution is about 
2-3 minutes in order to capture the details of the transition of the shock. 

The modeled fictitious ICME is launched on 15 February 2011 01:00:00 into 
a steady state background solar wind {Vsw ~ 430 km s""'^). It has a very wide 
angular extent within the ecliptic and is expanding to all sides. It is directed 
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toward the virtual in-situ spacecraft at the position of Earth and the ICME 
apex arrives 48 hours later on 17 February 2011 01:00:00. The ICME is faster 
than the ambient solar wind. It decelerates from 1100 km s~^ at 60 R© to 800 
km at 150 Rq, and the shock finally arrives with « 650 km s~^ at the in-situ 
spacecraft. 

Figure [T] shows the synthetic in-situ data extracted at the position of Earth. 
In the upper panel the three components of the magnetic field and the total 
magnetic field strength (black line) are plotted in cartesian coordinates with 
origin in Sun-center, x- and y-axes lie in the ecliptic plane and z-axis points 
to the ecliptic north pole to complete the right hand triad. The middle and the 
bottom panel show the speed in km s~^ and the proton number density in cm"'^, 
respectively. The red vertical line indicates the arrival of the shock preceding the 
ICME. If the shock is absent, the arrival time cannot be defined unambiguously. 
This could yield an additional uncertainty. 

The left part of Figure [2] shows the location of the observing remote space- 
craft (yellow dot) and in-situ spacecraft (red dot) as viewed from top at one 
time step of the ICME evolution. The black semicircle is the Thomson surface 



(TS: Vourlidas and Howard 2006 Howard and DeForest 2012) and the dashed 
yellow line is the LoS along the ICME front, i.e. the projection of the TS onto 
the ecliptic plane, from the view of the remote spacecraft for which the Jmap 
is produced. The corresponding Jmap is shown on the right side of Figure [2j 
The bright features are the integrated signal along the LoS with its intensity 
maximum at the TS. The method to produce synthetic Jmaps is described in 



Lugaz et al. (2009b). 



3. Analysis Methods 



To determine the kinematical properties of an ICME, we usually measure the 
apparent front, either in the direct white-light images or within Jmaps. Tracking 
the bright features in a Jmap has two advantages. First, the contrast between 
dense and less dense areas can easily be seen, which enables us to track ICMEs 
over larger elongations, hence, radial distances. Especially in HI2 images the 
brightness rapidly decreases and it is hardly feasible to identify the front at 
larger elongations. Second, due to the shape of the track that is directly visible 
in the Jmap, i.e. the evolution of elongation with time, it is possible to get a first 
estimate of the direction of the ICME (Sheeley et aZ.l|2008p. This is used in fitting 
methods as the Fixed-^ fitting (iKahler and Webb / 2007 ; Rouillard et aZ. 1 120081) 



or the Harmonic Mean fitting methods (Lugaz |2010(|Mostl et al. 2011 ) that can 
also be used to predict ICME arrival times. However, these fitting methods use 
the assumption of constant speed, which is very useful to get a rough estimate for 
the direction and the arrival time of the ICME. But to analyze variations in the 
kinematical behavior due to interactions with e.g. high speed streams, corotating 
interaction regions or preceding ICMEs 



[e.g. 



Gopalswamy et al., 2001 



Plunkett, and St. Cyr[ |2002[ |Lugaz, Vourlidas, and Roussevt |2009a 



Burlaga, 



Temmer 



et al. 2011 2012) it is required to apply a method, which allows the ICME 



speed to vary during propagation. 
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Figure 1. Synthetic in-situ data at the position of Earth. The upper panel shows the magnetic 
field with its three components (cartesian coordinate system) and the total magnetic field 
strength, the middle panel shows the speed and the bottom panel the proton number density. 
The vertical red line marks the shock arrival time at Earth that is used for further calculations. 
The X-axis gives the time in hours from launch of the CME at the Sun. 
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Figure 2. Left: top view of the modeled ICME for 60° separation between remote (yellow 
dot) and in-situ spacecraft (red dot). The black semicircle is the Thomson surface and the 
dashed yellow line is the LoS along the ICME front. Right: side view of the propagating ICME 
simulated as observed by STEREO-A/HI. The bright track shows the evolution in time of the 
apparent ICME front. The black cross marks the elongation of the in-situ spacecraft at shock 
arrival time. 
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Such an approach, which uses additional in-situ data to constrain measure- 
ments from white- hght data, is described in Rollett et al. ( 2012[ ). In this study 
we vaUdate the CHM method (using the constraint of arrival time) and derive 
the effective error in the ICME distance-time and speed-distance profiles. 

3.1. Constrained Harmonic Mean Method 



In the following we briefly outline the CHM method and its main features. 



for more details we refer the interested reader to Rollett et al. (2012). At this 



point we stress that the CHM method is not a forecasting method but an 
approach to analyze ICMEs after they are remotely observed and measured 
in-situ^ respectively. In this way we are able to derive a kinematical profile for 
the purpose of analyzing variations in speed. In order to convert the measured 
elongation angle [e(i)] into a unit of distance we use an approach assuming a 
circular ICME front, the so-called Harmonic Mean method (HM: Howard and 
Tappin 2009 Lugaz, Vourlidas, and Roussev 2009a I. HM was shown to be 



a quite good approximation for very wide events. This method assumes the 
observer to look along the tangent to the spherical ICME front and more or less 
disregards the concept of the intensity maximum on the TS. Equation ([l]) shows 
the conversion of e{t) into the radial distance of the ICME apex [i?HM(^)] for a 
given propagation direction [0]: 



RmA{t) 



(1) 



2do sin e(t) 
1 + sin(e(i) + 4>) ' 

where do is the Sun-observer distance. For the constraint with the in-situ arrival 
time we need the distance of the ICME front in the direction of the in-situ 
spacecraft. This is given by 



RmAiu W = ^HM,ta (</>) cos S, (2) 

where i?HMi,ta(0) is the distance of the ICME front in the direction of the in-situ 
spacecraft at arrival time [ta], Rnu.t^ is the distance of the ICME apex at arrival 
time and S is the angle between the line of the in-situ spacecraft to the Sun [di] , 
and Rauit)- To determine the propagation direction we now calculate 



AdHMi(0) = ^HMi,t,(0) - di, 



(3) 



for (j) e [0°,180°]. The direction providing the minimum value for AduMiW 
is the resulting propagation direction (j) that is now used to convert e{t) into 
distance by using Equation ([T]). 

In other words: the input parameter (j), and therefore also the shape of the 
distance-time profile, is changed until the best accordance with the in-situ arrival 
time is found. This ensures that the HI observations are consistent with the in- 
situ data, which serve as a boundary condition for the kinematical profiles. In 
addition, the input paramter ip, used for the best agreement, is the calculated 



propagation direction of the ICME apex (Mostl et al. 2009b 2010 Rollett et al 
Mil). 
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3.2. Speed Calculation 



When the ICME speed profile is derived from the calculated distance its shape 
strongly depends on the uncertainties in the distance profile. Therefore it is 
necessary to fit the distance profile and use this fit for further calculations. In 
Gallagher, Lawrence, and Dennis (2003) and'Bein et al. (2011) fitting functions 



and methods are described for coronagraphic FoVs close to the Sun. These func- 
tions take into account the strong acceleration and deceleration at the beginning 
of the CME evolution (Zhang et al. 20011. Since we analyze data from HI we 
are mainly interested in the propagation phase, during which the ICME adjusts 
to the solar wind speed, and developed a function, which takes into account the 
different parts of the distance-profile for the HI FoV: 



F{x)^ Aihi{x)+A2x + A'i^/^ + Ai, (4) 

where the first term reproduces the strong deceleration at the beginning of the 
HIl FoV, the second term describes the roughly constant speed during the propa- 
gation phase. The third term takes into account the slowly decreasing/increasing 
speed and is needed to avoid an overestimation of the arrival speed at 1 AU, and 
the constant is the y-intercept. Figure [3] shows the distance-time fits and the 
fitting parameters for all three separations. The time derivatives of the distance- 
time fits are over plotted on the speed-time profiles, derived from the distance 
measurements converted from the elongation-time profiles. The error bars show 
the standard deviations calculated from measuring each track five times. 



3.3. Error Analysis 



Uncertainties in the analysis of ICME kinematics and directions are due to 
various reasons. By using a modeled ICME it is possible to provide an error 
estimation for the most significant contributions. 

If the apparent front of an ICME is measured from white-light images the 
intensity range is varying, i.e. if a track is measured by more than one person 
there will be a certain spread of measurements. To quantify this error we de- 
termine an error range for the measurements by performing two surveys — one 
very high (outermost) and one very low (innermost) track (see Section 4.1 ). The 
mean differences between these two tracks lie between ±0.3° (HIl) and ±1.35° 
elongation (HI2) for all three separation cases. 

To estimate the error due to visually tracking of the elongation-time profile, 
the Jmap track is measured five times. Out of this sample we calculate the mean 
values and standard deviations for every Jmap track. The elongation errors for 
all three separation cases lie in a range of about ±0.1° (HIl) and ±0.4° (HI2). 

Table [T] lists the mean values for the three cases of different separation angles 
between in-situ and remote spacecraft and both HI instruments, for the two de- 
rived errors in elongation, i.e. the intensity range, where measurements could be 
taken within the Jmap and the measurement error itself. For further calculations 
only the measurement error is taken into account. Table [2] gives the errors for 
the resulting distance in Rq and the speed in km s~^. 
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Figure 3. Distance-time profiles with over plotted fits for 30°, 60°, and 90° separation of 
remote and in-situ spacecraft. In eacfi case tlie upper panel shows the distance measurements 
converted from the elongation-time profiles (triangles), the over plotted fits and the fitting 
parameters in Rq . The lower panels show the time derivatives from the distance-time profiles 
and the over plotted time derivatives of the distance-time fits. The red crosses mark the in-situ 
measured arrival time and speed, respectively. 
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We only consider errors due to the measurements itself. Uncertainties bea- 
cause of the assumed circular shape or the angular width of the ICME front are 
not taken into account. Since HM assumes a width of 180° it may not be the 
right approach for narrow events. In this case the Ftjj method should provide 
more reliable results. 

Table 1. Summary of calculated errors denoted in elongation for each separation 
case and each heliospheric instrument: HIl and HI2. The first row lists the standard 
deviation of the innermost and the outermost Jmap track, i.e. the possible Jmap 
tracking area. The second row shows the measurement error derived out of the five 
tracks in elongation. 

30° separation 60° separation 90° separation 
HIl HI2 HIl HI2 HIl HI2 



Jmap tracking area ±0.25° ±1.35° ±0.3° ±1° ±0.3° ±0.7° 

mean standard deviation ±0.1° ±0.4° ±0.5° ±1.9° ±0.5° ±1.4° 



Table 2. Resulting errors for the ICME kinematics. The first row lists the 
distance error in Rq based on the mean standard deviations for both in- 
struments. The second row shows the speed error derived from the standard 
deviations. 





30° 


separation 


60° 


separation 


90° 


separation 




HIl 


HI2 


HIl 


HI2 


HIl 


HI2 


distance error [Rq] 


0.4 


1.7 


0.4 


1.6 


0.4 


1.3 


speed error [km s^-*^] 


31 


128 


34 


102 


27 


81 



4. Results 

In this section wc provide the calculated directions and the kinematical pro- 
files derived from the measurements within the synthetic Jmaps by using the 
CHM method. For validation, these results are compared to the "true" apex 
kinematics, measured from the top view of the modeled ICME. 

4.1. Kinematics and Direction 

Due to the fact that the ICME is directed toward the virtual in-situ spacecraft, 
that is located at the position of Earth, the real ICME direction is EO/WO. For 
all three separation cases the Jmaps with the outermost and innermost tracks 
are shown in Figure [4] The red line is the innermost and the blue line is the 
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outermost edge of the ICME front. The black cross marks the elongation of the 
in-situ spacecraft at shock arrival time. Because of the large angular width of 
the ICME, LoS effects make the ICME to appear earlier at 1 AU than actually 
measured in-situ. 

The resulting directions obtained for the mean of the innermost and the 
outermost Jmap tracks are shown in Figure [Sj The yellow circle is the assumed 
circular shape of the HM approximation and the yellow arrow marks the direction 
of the ICME apex. The dashed semicircle is the TS. For all three separation cases 
the resulting directions are consistent within a few degrees (30° sep: = W7; 60° 
sep: 4> = W12; 90° sep: = W15; true dh.: EO/WO). Because of the corotating 
interaction region (CIR), heading in front of the ICME, the western flank is 
denser, which may be the reason for a result toward the western direction. The 
resulting propagation angles for the innermost tracks are 2-4° larger than the 
directions calculated from the outermost tracks (see Table |3| . To demonstrate 
how well the CHM method works, we calculate the directions also for the con- 
strained F(j) method using the constraint of the in-situ data point of the arrival 



time of the ICME (see [RoUett et al. 2012). Fcj) assumes a single point moving 



radially away from the Sun and is thus more dependent on the separation of 
both spacecraft and/or the remote spacecraft to the ICME apex than HM. 

Table 3. Resulting ICME propagation directions for the three 
separation cases. Constrained Fixed-(/) (CF0) and constrained 
Harmonic Mean (CHM) methods are used to calculate the di- 
rections for the innermost (inn.) as well as for the outermost 
(out.) track in the Jmap. 

30° separation 60° separation 90° separation 



CHM 
CF(f> 



inn. 


out. 


inn. 


out. 


inn. 


out. 


W6 


W8 


WIO 


W14 


W15 


W16 


W13 


W17 


W22 


W28 


W32 


W25 



In Figure |6] the converted Jmap measurements (solid lines) are compared to 
the real ICME apex distance-time profile (dashed line) . The top panel shows the 
30° separation case which overestimates the distance especially for the HIl FoV. 
The best match of measurements and ICME apex is given for the innermost edge 
(red line) in the Jmap. The middle panel shows the distance-time profile for 60° 
separation revealing that both tracks are in a relatively good agreement with the 
ICME apex compared to the other separation cases. The bottom panel shows 
the 90° separation case. In contrast to the 30° separation case, the distance is 
underestimated most of the time. For 60° and 90° separation the measurements 
taken at the outermost edge of the Jmap track reveal the best agreement to the 
apex measurements. 

Figure [7] displays the same approach as Figure |6] but for the speed-distance 
profile. The ICME apex enters the HIl FoV with a speed of 1100 km s~^. 
It slowly decreases to a mean speed of w 750 km s^^. The "true" speed profile 
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Figure 4. Jmaps for the three separations with the measured tracks of the ICME front over 
plotted. The red line is the innermost track and the blue line is the outermost track. The black 
crosses mark the elongation of the in-situ spacecraft at shock arrival time. 
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Figure 5. Top view of the modeled ICME. The denser region (red part) is located at the 
western flank of the ICME. The spacecraft positions (yellow and red dots) and resulting 
directions from the CHM method (yellow arrow) for all three separations (top: 30°, middle: 
60°, bottom: 90°) are over plotted. The black semicircle is the TS and the yellow dashed line 
is the tangent to the circle assumed by the Harmonic Mean approximation. The blue cross 
marks the ICME apex that was measured from top view. 
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Figure 6. Distance— time profiles for 30°, 60°, and 90° separation from the in-situ spacecraft. 
The dashed Une shows the top view measurement of the apex, the solid lines are the conversions 
from the innermost (red) and the outermost (blue) Jmap measurements., the shaded areas are 
the error ranges. The crosses indicate the arrival time at the in-situ spacecraft. 
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of the ICME apex is derived from the top view images. With growing distance 
from the Smi its standard deviation decreases from 148 km s^^ to about 62 
km s~^ because of the increasing resolution of the simulation. The speeds de- 
rived from the innermost Jmap tracks for all three separation cases show a 
stronger deceleration than the speed profiles from the outermost edges. The 
30° separation case reveals the best agreement, while the 60° as well as the 90° 
observing angle in the early propagation phase under- and in the late propagation 
phase overestimate the ICME speed. With growing separation between remote 
spacecraft and ICME apex the overestimation of the arrival speed increases. 
While in the 30° separation case the arrival speed of the outermost edge is 
quite close to the apex speed, the innermost track underestimates the arrival 
speed by about 100 km s~^. From both other viewing angles the arrival speed 
is overestimated on average by about 75 km s~^ (60° separation) and 125 km 
s^-^ (90° separation), respectively. As well as for the distance-profile, measuring 
the Jmap track at the innermost edge in the 30° and 60° separation cases gives 
more consistent results compared to the apex. 



5. Discussion and Conclusions 

One aim of this study was to validate the constrained Harmonic Mean method 



(CHM) introduced in Rollett et al. (2012). This method assumes a spherical 



ICME front always attached to the Sun and the observer to look along the 
tangent of the front. The CHM method uses single-spacecraft heliospheric imager 
(HI) observations and in-situ data at 1 AU. By constraining the kinematics 
with the in-situ arrival time it is possible to derive the propagation direction for 
the ICME apex. In addition we can provide distance-time and speed-distance 
profiles for the entire propagation phase. The second purpose was to provide a 
reliable error estimate for the kinematics of the ICME when applying the CHM 
method to Jmap measurements. 

In order to test the CHM method we used a numerical simulation of an ICME 
and analyzed synthetic Jmaps from three different vantage points (30°, 60°, 
and 90° separation of remote spacecraft and in-situ spacecraft). The resulting 
directions and kinematics for each of the three separation angles were compared 
to the "true" kinematics of the ICME apex, i.e. free from effects due to Thomson 
scattering or conversion methods, derived from the top view of the modeled 
ICME. 

The directions obtained seem to be not very sensitive to the separation angle 
between remote and in-situ spacecraft (30° sep: (j) = W7; 60° sep: (j) = W12; 
90° sep: (/> = W15; true dir.: EO/WO). The density of the ICME is higher in 
its western flank as a result of a CIR ahead of the ICME. This fact could 
have influenced our results for the propagation direction and shifted them to 
a western direction. Another possible reason for that could be that the CHM 
method assumes a smaller curvature radius of the ICME front than shown by the 
simulated ICME. It still has to be tested, how the results of the CHM method 
depend on the ecliptic extend, i.e. the width of the ICME. 
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Figure 7. Speed-distance profiles for 30°, 60°, and 90° separation from the in-situ spacecraft. 
The dashed line shows the time derivatives from the top view distance measurement of the 
apex, the solid lines are the time derivatives from the conversions from the innermost (red) 
and the outermost (blue) Jmap measurements. The crosses indicate the arrival speed at the 
in-situ spacecraft. 
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Future studies should prove the behavior of the method when ICMEs with 
different widths are investigated. In Lugaz, Roussev, and Gombosi (2011) the 
results of the F(/) and HM fitting methods, which assume constant speed, were 
tested by analyzing a simulated ICME from different viewing angles. Consistent 
with our study, a stronger dependance of the observing angle for F(j) than for 
HM was obtained when calculating the propagation direction. 

The ICME kinematics were compared to the actual distance-time and speed- 
distance profiles. The best agreement for the distance is obtained for the 60° 
separation case. The 30° separation case overestimates the distance most of 
the time while in the 90° separation case the distance is too small compared 
to the ICME apex derived from the top view. For the speed-distance profiles 
the best accordance is revealed by the 30° separation case. Both other settings 
first under- and then overestimate the speed. The differences between the true 
and the calculated arrival speed increases with growing separation (30° sep: 
AVarr « "50 km s-i, 60° sep: AKit « +75 km s-\ 90° sep: AVarr « +125 
km s"""^). In general, conversion and fitting methods making special geometrical 
assumptions on the shape of the ICME front, e.g. Fcj), HM or SSE methods, 
overestimate the arrival speed with increasing observation angle. 

The error estimation was done by taking into account two different error 
sources. The intensity area in the Jmap, where the apparent ICME front could be 
measured, varies between ±0.3° (HIl) and ±1.35° elongation (HI2). This should 
reflect the range of tracks if more than one person measures the front. This is in 
agreement with the results found by Williams et al. { 2009 1 who quantified the 
elongation range of chosen tracks measured within a Jmap with less than ±1° for 
HIl and ±2° for HI2. We find that in general, measuring the Jmap track at the 
innermost edge gives more consistent results compared to the apex kinematics. 
The second error is due to the visual measuring of the front and is obtained by 
calculating the standard deviations out of five measurements taken that lie in a 
range of ±0.1° (HIl) and ±0.4° (HI2). The errors within the speed-profiles lie 
in a range of 27-34 km s^^ (HI) and 81-128 km s^^ (HI2). 

The two STEREO spacecraft already are on the backside of the Sun, and thus 
a forecast for Earth directed ICMEs using the heliospheric imagers is no longer 
possible. To improve existing space weather forecast an L4/L5 mission carrying 
heliospheric imagers is a valuable option. The side view of the entire distance 
from the Sun up to 1 AU is a big advantage upon coronagraphic observations 
from Earth view. 
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